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We study the two dimensional quantum Heisenberg antiferromagnet on the square lattice with easy- 
axis exchange anisotropy. By the semiclassical method called pure-quantum self-consistent harmonic 
approximation we analyse several thermodynamic quantities and investigate the existence of a finite 
temperature transition, possibly describing the low-temperature critical behaviour experimentally 
observed in many layered real compounds. 

We find that an Ising-like transition characterizes the model even when the anisotropy is of the 
order of 10~'^J (J being the intra- layer exchange integral), as in most experimental situations. On 
the other hand, typical features of the isotropic Heisenberg model are observed for both values of 
anisotropy considered, one in the guasi-isotropic limit and the other in a more markedly easy-axis 
region. 

The good agreement found between our theoretical results and the experimental data relative to 
the real compound Rb2MnF4 shows that the insertion of the easy-axis exchange anisotropy, with 
quantum effects properly taken into account, provides a quantitative description and explanation of 
the experimental data, thus allowing to recognize in such anisotropy the main agent for the observed 
onset of finite temperature long-range order. 



I. INTRODUCTION 

In recent years the increasing interest in the physics 
of low-dimensional magnets has led to a deep analysis of 
the thermodynamic behaviour of two-dimensional quan- 
tum antiferromagnetsci and, in particular, of the isotropic 
quantum HeiS|enberg antiferromagnet on the square lat- 
tice (QHAF)0l3. Such model has in fact been widely used 
to describe the magnetic properties of many quasi two- 
dimensional real compounds, from the S = 1/2 cuprate 
La2Cu04QEl to the S = 5/2 compound Rb2MnF4, re- 
cently studied by Lee et al.cl. _ |_. 

Different theoretical methodaOilj have been used to ex- 
amine the rich reservoir of experimental data and the pic- 
ture of the subject is now well focused, albeit with pc|raej 
shaded parts. In particular open questions still existE5~t£l 
on the low-temperature region, where the spin correlation 
length becomes of the order of 10^ lattice spacings and 
the real magnets are seen to develop macroscopic areas 
of correlated spins. 

The theoretical debate on the low-temperature regime 
has been mainly dedicated to the isotropic QHAF, but 
real compounds are not actually well described by such 
isotropic model when temperature is lowered: the exper- 
imental evidence of a finite-temperatuBei transition, op- 
posite to the Mermin- Wagner theoremt^l assertion that 
such a transition cannot occur in the two-dimensional 
isotropic QHAF, suggests that three-dimensional corre- 
lations and possible anisotropy effects, as well as a com- 
bination of both, must be considered. 

The magnetic structure of layered real compounds 
is such that the exchange integral J for neighbouring 



spins belonging to the same plane is orders of magni- 
tude If^pi than that for neighbouring spins on different 
planeg23~E3, hereafter called J'; one would hence naively 
expect the magnetic properties to be those of an effec- 
tive two-dimensional magnet down to temperatures of 
the order of J', until the transition towards an ordered 
three-dimensional phase should take place. However, the 
experimentally observed transition occurs at a critical 
temperature of the order of J, signalling the transition 
itself to be driven by the intra-layer exchange interac- 
tion; details of such interaction, such as possible easy- 
axis or easy-plane anisotropics, are hence fundamental in 
the analysis of the critical behaviour. 

Several works (see Ref. ^ for a review) have shown 
that, besides the superexchange interaction, there ex- 
ist further interaction mechanisms whose effects may be 
taken into account by inserting proper anisotropy terms 
in the magnetic Hamiltoiiian; in particular, JJae transi- 
tion obsecyed in KaNiEjtB (S = 1), Rb2FeF4B (5 = 2), 
K2MnF4N, Rb2MnF4E3 {S = 5/2), and others, is seen 
to be possibly due to an easy-axis anisotropy. Such 
anisotropy has often been described in the literature via 
an external staggered magnetic field: This choice, de- 
spite allowing a qualitative description of the experimen- 
tal data, lacks the fundamental property of describing a 
genuine phase transition, as the field explicitly breaks the 
symmetry and makes the model ordered at all tempera- 
tures. 

In order to produce a second order transition known to 
be due to an easy-axis anisotropy, it is actually appropri- 
ate to insert such anisotropy in the form of an exchange 
one, thus preserving the symmetry under inversion along 
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the easy-axis, whose spontaneous breaking manifests it- 
self in the transition: subject of this paper is the study 
of the thermodynamic properties of the resulting model, 
hereafter called Easy- Axis Quantum Heisenberg Antifer- 
romagnet (EA-QHAF). 

In Sec-H we define the model and discuss its general 
properties, posing the problems we want to address. The 
method used is briefly described in Sec.p|, where results 
for several thermodynamic quantities and different values 
of the anisotropy are shown and commented. In Sec. IV 
we compare our results with the available experimen- 
tal data for the staggered magnetization and suscepti- 
bility, and for the coHrglation length, of the S = 5/2 real 
compound Rb2MnF4la'E2l, while conclusions are drawn in 
Sec.0 



II. MODEL AND METHOD 



harmonic approximation (PQSCHA) already suc- 

cesfully applied to many magnetic systems in one and two 
dimensions. The method reduces quantum expressions 
for statistical averages to effective classical-like ones, con- 
taining temperature and spin-dependent renormalization 
parameters. The thermodynamics of the effective model 
can then be studied by means of classical techniques, like 
Monte Carlo simulations. 

The PQSCHA is known to be particularly suitable for 
anisotropic (easy-plane or easy-axis) models, but it also 
gives very good results in the isotropic QHAF so that we 
can confidently use it to investigate the possible crossover 
between a Heisenberg-like behaviour at high temperature 
and an Ising-like one near T^, such crossover is detected 
in the experimental data to such an extent that, well 
above the transition, real compounds can be satisfacto- 
rily described by the isotropic QHAF . 



The EA-QHAF on the square lattice is described by 
the Hamiltonian 



III. RESULTS 
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where i = (ii, 12) runs over the sites of a square lattice, d 
connects each site to its four nearest neighbours, J > 
is the antiferromagnetic exchange integral and fj, is the 
anisotropy parameter (0 < < 1 for easy-axis models). 
The spin operators S" (a — x,y, z) are such that |Sp — 
S{S + 1) and obey [S*?, S?] = ie^^p^SviS^^ . 

When /Lt = 1 the model looses its easy-axis character 
and reduces to the isotropic QHAF. 

The /X = case will be hereafter called Ising limit, not 
to be confused with the genuine Ising model, reproduced 
by Eq.(|l|) with /i = and S — 1/2. Despite being a 
very particular case of Eq. (|l|) , the two dimensional Ising 
model on the square lattice is a fundamental point of 
reference for the study of the thermodynamic properties 
of the EA-QHAF. 

In particular, a renormalization-group analysis of the 
classical counterpart of the model Eq.(|l|)EZE3 foresees 
the occurrence of an Ising-like transition at finite tem- 
perature for any value of /i, no matter how near 
to the isotropic value /i = 1; this analysis is sup- 
ported by seyetal works based on classical Monte Carlo 
simulation£D"E3. 

In the quantum case, however, no information is given 
about the value of the critical temperature Tc{^,S) as 
a function of the anisotropy rS-pd of the spin, save the 
fact that Tc(0,l/2) = 0.567 ^1 and rc(l,S') = 0. We 
do not hence know whether or not the small anisotropy 
(l/i — 1| ~ 10~^) observed in real compounds can be re- 
sponsible of transitions occurring at critical temperatures 
of the order of J, given also the fact that we expect quan- 
tum fluctuations to lower the critical temperature with 
respect to the classical case. 

We have developed a quantitative analysis of several 
thermodynamic properties of the model, by means of the 
semiclassical method called pure-quantum self-consistent 



The main output of the PQSCHA is the effective 
Hamiltonian appearing in all statistical averages, whose 
expression for the model described by Eq.(^ reads 

i,d 



(2) 



where s — (s^, s^, s^) is a classical unit vector, 5 = 5-1- 
1/2, and t = T / JS"^ is the reduced temperature hereafter 
used. The appearance of the minus sign in front of the 
effective Hamiltonian is due to the fact that in Eq. (||) , as 
in all the classical-like expressions reported below, spins 
belonging to one of the two sublattices have been flipped, 
being this an innocuous operation at a classical level. 
The renormalization parameters are 
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where the coefficients 
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are self-consistently determined by solving Eqs.(||) and 
(0), with 



(5) 



and 
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k being the wave vector in the first Brillouin zone. 
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The temperature and spin dependent uniform term 

/ sinh /k ' 

111 

k 



g{t,s) = 2js'^{i-epl) + js^tJ2^n 



does not enter the expressions for the statistical averages, 
but contributes to the free energy and to the related ther- 
modynamic quantities. 

The renormalization coefficients and measure 
the pure-quantum fluctuations parallel and perpendicu- 
lar to the easy-axis, respectively, of one spin with respect 
to its nearest neighbours, and vanish in both the high- 
temperature and the classical S oo limit, ft is worth- 
while noticing that the fluctuations along the easy-axis 
only enters the renormalization of the z-component of 
the exchange interaction. 

Defining the effective exchange integral and effective 
anisotropy 



JoS = J 



el 

Meff = M 



Eq.(0) can be written in the form 
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Heft = ~^JcsS^ J2 [^-ff i-'i^i+d + ^f^f+d) + Sf+d] 
i.d 



5it, S) 



(6) 



to make evident that the PQSCHA leads to a classical- 
like effective model of the same form of the original 
quantum one, whose thermodynamic properties can be 
studied by classical numerical technique, properly taking 
into account the spin and temperature dependence of the 
renormalized parameters Jeff and ficS- As 0y < 6*^ < f , 
it is Jeff < J and f > ficS > so that quantum effects 
are seen to cause the weakening of both the exchange in- 
teraction and the easy-axis anisotropy; in the isotropic 
limit — T>± and ficS = ^ = f . 

The PQSCHA expression for the statistical average of 
a quantum operator O is 



:o) = 
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d"s O. 
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where f3 = T-\ Z^g = J d"se-^^<=ff; Oeff is deter- 
mined by the same procedure leading to HeS and contains 
temperature and spin dependent renormalizations; after 
Eq.([7|) most thermodynamic quantities may be written 
in a particularly suggestive form in terms of classical-like 
statistical averages defined by the effective Hamiltonian 
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as far as the evaluation of ( • ■ • )eff is concerned, one must 
keep in mind that, because of the temperature depen- 
dence of the effective anisotropy appearing in Tieff, each 
point in temperature corresponds to a different effective 
model. This means that, if the classical Monte Carlo 
technique is used, as done in this work, the simulated 



model changes with temperature, so that, at variance 
with the isotropic case, no existing classical data can be 
used and a complete series of ad hoc simulations must 
be carried on. Nevertheless, the computational effort re- 
quired is still that of a classical simulation, as the evalu- 
ation of the renormalized parameters is a matter of a few 
seconds on a standard PC. 

The apphcation of the PQSCHA to the EA-QHAF 
leads to the following results for the thermodynamic 
quantities we have considered: 



• Internal energy 
u={n)/{NJS^): 



u ' 
where 



NJS- 



■{HcE)os~g{t,s) + Ht,s) 



(8) 
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is a negative zero-point quantum correction term 
due to the anisotropy (it vanishes in the isotropic 
limit). 



• Staggered magnetization 
m EE X;i(-l)''+'M'S'f)/iV5: 

m-0i(s?)eff + M(t, S) 



(9) 



where the magnetization renormalization is seen to 
be due to an effective spin reduction {9\ < 1)) 
and to the appearance of the negative term A4 = 
— {Ol — 6'p/(l — /i^) which vanishes in the isotropic 
limit and is finite for fi — 0. 

Staggered correlation function above Te G{r) = 
(-f)''i+'~^(Si-Si+r)/52 withr = (ri,r2) any vector 
of the square lattice and r = |r|: 



G{r) = 0^{Si ■ Si+r)eff , 



(10) 



where 0^ ^l-V,,/2 and 



- cos(k-r)) £k 



is a further site-dependent renormalization coeffi- 
cient. For increasing r, the coefficient 2?r rapidly 
converges to a uniform term, so that the asymp- 
totic (r — > oo) behaviour of G(r) is actually deter- 
mined by that of the effective classical-like correla- 
tion function (Si • Si-|_r)eff- 



• Staggered static susceptibility above Tc 
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• Correlation length above Tc We have determined 
the correlation length ^ by fitting G(r) with 
the expression proposed by Serena, Garcia and 
LevanyukE^ 



G(r) oc 



1 



^4 (^/^)^ + (^/^). 



which interpolates the two asymptot 
r — > oo and r — *■ of the Ising model 



(12) 



behaviours 



In what follows, we show our results as obtained com- 
bining the above PQSCHA expressions with the numer- 
ical output of the classical Monte Carlo simulations we 
have performed to evaluate the efTective statistical aver- 
ages ( ■ ■ • )off - At variance with the isotropic case, where 
results for different values of the spin are obtained by 
the same series of classical simulations, we now have to 
fix the value of the spin in order to determine, for a 
given value of /i, the corresponding /ioff to be used in 
the simulation. We have hence concentrated ourselves 
on the case /i = 0.9942 and S — 5/2, being these the 
anisotropy and spin values corresponding to the real com- 
pound Rb2MnF4; the more anisotropic case fi = 0.7 and 
5 = 5/2 has also been considered, because of its expect- 
edly more marked Ising-like features. 
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FIG. 1. Internal energy vs t for S — 5/2 and /i — 0.7 
(dashed), 0.9942 (full), 1 (dash-dotted); the curve relative to 
the Ising model {S = 1/2 and fi — 0) is shown in the inset. 

In Fig.|l| and ^ we show the internal energy and specific 
heat, versus temperaturp- for both values of compared 
with the isotropic casetHi ^ = 1. As the value of the 
critical temperature increases for larger anisotropy, but 
diminishes when smaller spin values are considered, we 
avoid the confusing direct comparison between our curves 
{fi > 0,S = 5/2) and the correponding quantities for the 
Ising model {fi = 0,5 = 1/2), by showing the latter as 
Insets. 

The Ising character of the /i = 0.7 model is evidenced 
by the pronounced peak in the specific heat, correspond- 
ing to a qualitatively different temperature dependence of 
the internal energy with respect to that of the isotropic 
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0.0 0.4 0.8 1.2 

FIG. 2. Specific heat vs t (lines and inset as in Fig.l 



model. Although such difference is almost not percep- 
tible when the fi = 0.9942 curve is considered, a clear 
peak in the specific heat is still present, testifying to a 
persistence of the Ising-like behaviour even in this quasi- 
isotropic model. 

The same conclusion is drawn when the staggered 
magnetization is considered: in FigJ^ we see that for 
fi = 0.9942 there exist a wide temperature range where 
the system is ordered, with a critical temperature tc that, 
despite being lower than the one relative to the = 0.7 
case, is still of the order of J. In the inset we show the 
magnetization curves normalized to their saturation val- 
ues, as functions of t/t^, together with that of the Ising 
model: It is evident that an increased anisotropy causes 
a sharpening of the way the magnetization vanishes. 
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FIG. 3. Staggered magnetization vs t for S = 5/2 and 
fj, = 0.7 (dashed), 0.9942 (full); the dotted lines are the data 
fits used to extract the critical temperature values. In the 
inset the same curves, normalized to the saturation value, 
are shown versus t/tc, together with the magnetization of the 
Ising model (dotted). 

The critical temperatures hereafter used are tc = 0.785 
for fi ~ 0.7 and tc ~ 0.575 for ^ = 0.9942; these values 
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have been determined by locating at best the correlation 
length divergence and consistently coincide with those 
emerging from the analysis of the critical behaviour of 
other quantities. For instance the critical temperatures 
determined by fitting the magnetization curves (see the 
dotted lines in Fig.y) are t'^ = 0.787 for ^ = 0.7 and 
t'^ = 0.576 for ^l = 0.9942. 

It is to be noticed that even below ic, where the finite 
value of the magnetization would suggest a complete pre- 
dominance of the Ising character, the system does actu- 
ally display features which are distinctive of the isotropic 
model. In particular, the specific heat for both values of /i 
shows, after an exponential start typical of a gapped dis- 
pertion relation, a change in the curvature and an almost 
linear temperature dependence, due to the excitation of 
long-wavelength low-energy modes of the same type of 
those characterizing the isotropic model. It is just in 
the vicinity of tc that a new change in the curvature an- 
nounces the forthcoming transition. 
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FIG. 4. Correlation functions for 5" = 5/2, = 0.9942 and 
t = 0.6, 0.7, 0.85 (from the top)). 

In Fig.^ we show the correlation function G{r) as a 
function of r in the quasi-isotropic case ji = 0.9942 
and for three different temperatures. The fit of our 
data with the Serena-Garcia-Levanyuk function Eq.(|l^) 
is very good in the whole temperature range examined, 
and hence leads to a clean evaluation of the correlation 
length. 

In Fig.H we show the correlatioiiJength and also report 
the curve for the isotropic modelo: we notice that the 
/i = 0.9942 curve lays on the isotropic one up to cor- 
relation lengths of the order of 20 lattice spacings (i.e. 
t ~ l.OStc), while for /i = 0.7 a deviation is evident 
already for ^ ~ 2 (i.e. t ~ 1.3ic)- This means that, 
in the former case, there is a temperature region where 
the model significantly behaves like the isotropic model, 
as far as the correlation length is concerned, and it is 
hence meaningful to introduce the idea of a crossover 
from a Heisenberg- to an Ising-like regimeQ; on the con- 
trary, when /i = 0.7 the Ising character of the model 
is manifest already when correlations over a few lattice 
spacings develop. 
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FIG. 5. Correlation length vs t for S — 5/2 and fi — 0.7 
(dashed), 0.9942 (full), 1 (dash-dotted). 



IV. COMPARISON WITH EXPERIMENTAL 
DATA 

The results shown in previous section qualitatively 
explain the mechanism possibly underlying the finite- 
temperature ordering experimentally observed in many 
quasi-two dimensional real magnets. 

As for a more precise quantitative analysis, we have 
concentrated ourselves on the 5 = 5/2 magnet Rb2MnF4: 
reason for this fihoice is the availability of recent neutron 
scattering datafl relative to such compound and the fact 
that, because of its cristallographic structure, Rb2MnF4 
is known to behave as a two-dimensiO|Ua.l .magnet both 
above and below the observed transitionE3E3. This means 
that the critical behaviour is not contaminated by the on- 
set of three-dimensional order and a clean characteriza- 
tion of the transition is possible, as well as a meaningful 
comparison with the experimental data for the magneti- 
zation below Tc. 

The model parameters Js — 7.&2,± 0.09 K and /is = 
0.9953 available in the literatureE3 for the compound 
Rb2MnF4 are obtained by fitting the extrapolated T ^ 
experimental data for the spin- wave frequencies with the 
expression 



(13) 



this means that Jg and fig are renormalized by the zero- 
point quantum fluctuations and are not the bare values 
to be inserted in Eqjl). These have hence been deter- 
mined equating Eq.(|l3) with the zero-temperature dis- 
persion relation relative to the EA-QHAF as given by 
the PQSCHA 



Wk = 4JS'/z6li(0)i 



/ ^1(0) 

M^^i(o) 



(14) 



where ^|(0) and 0^(0) are the renormalization parame- 
ters defined in Eq.(P) evaluated at t = 0. 
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The resulting equation for fi 

must be self-consistently solved, as both 9± and 0|| de- 
pend on /i, and gives /i — 0.9942. 

Once fi is determined, the equation for the exchange 
integral 

S Js 

is straightforwardly solved and gives J = 7.42 K; 
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FIG. 6. Staggered magnetization vs t for Rb2MnF4, nor- 
malized to the saturation value mo: our results (full line) are 
compared with experimental data (full circles) from Ref. |2^; 
the dashed line is the interpolating curve therein proposed. 

In Fig.^ our results for the staggered magnetization 
are shown together with the experimental data from Ref. 
^ and the interpolating curve there proposed. Besides 
the overall agreement in the whole ordered phase, it is to 
be noticed that our prediction for the value of the criti- 
cal temperature perfectly coincides with the one deriving 
from the experimental analysis, which gives Tc = 38.4 K 
(i.e. tc — 0.575); such a large value of Tc with respect to 
the exchange integral J — 7.42 K, should not surprise, as 
the squared of the spin value has been actually extracted 
from the latter. 

In order to obtain the best quantitative description of 
the EA-QHAF in the paramagnetic phase, and given the 
small anisotropy of Rb2MnF4, the lowest-temperature 
data for the staggered susceptibility and correlation 
length, shown in Figs.0 and are determined by the 
PQSCHA version introduced in Ref. ^ and there shown 
to be the most appropriate to study the isotropic model. 
The difference between such version and the original 
one described in Sec.|| consists in the appearance of 
the renormalization parameters ^ , instead of ^| , in 

Eqs.(|), with = 6ijj j^-X»Jj|j_/2, Vf^^ being the renor- 
malization coefficients determined by the classical self- 
consistent harmonic approximation. 
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0.4 0.8 1.2 

t 

FIG. 7. Staggered susceptibility vs t for S = 5/2, 
— 0.9942 (full) and fi — 1 (dashed-dotted); symbols are 

neutron scattering data on Rb2MnF4 from Ref. ^. A zoom of 

the critical region is shown in the inset. 

The agreement between our results and the experimen- 
tal data is indeed noticeable, given also the fact that no 
best-fit procedure has been used. In addition, our results 
represent a clear improvement with respect to those com- 
ing from a mean-field treatment of the anisotropy, pro- 
posed by Keimer et alB, which enables one to derive the 
correlation length of the anisotropic model directly from 
the data of the isotropic one, as done in Ref. ^pterting 
from the PQSCHA results for the isotropic modeltHi itself. 
In particular, it is evident from FigJ| that the mean-field 
approach, apart from accounting for the existence of the 
phase transition, leads to an overestimate of the critical 
temperature, while the model with exchange anisotropy 
gives a very accurate estimate of tc- 

V. CONCLUSIONS 

In this work we have studied the easy-axis quan- 
tum Heisenberg antiferromagnet on the square lattice by 
means of the pure-quantum self-consistent harmonic ap- 
proximation; expressions for several quantum statistical 
averages have been determined for the general model, 
with any value of the spin and of the anisotropy. The 
numerical work, consisting in classical Monte Carlo sim- 
ulations on a properly renormalized model, have been 
concentrated on the S = 5/2, /i — 0.7 and S = 5/2, 
fi = 0.9942 cases, the latter corresponding to the real 
compound Rb2MnF4. 

We have shown that a finite temperature transition is 
present in both cases and that such transition is cleary 
of an Ising type; the value of the corresponding criti- 
cal temperature have been determined by the analysis 
of the temperature dependence of the correlation length 
and has been always found to perfectly agree with that 
extracted by the analysis of other thermodynamic quan- 
tities. 

Despite the essential presence of the transition, the 
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FIG. 8. Correlation length vs t; lines, symbols and inset as 
in Figj^ apart from the dashed curve, representing the result 
of a mean-field approach to the anisotropy proposed in Ref. 
^ and the triangles, which are Quantum Monte Carlo data 
for the isotropic model, from Ref. |l6[ 



EA-QHAF does also display, both below and above tc, 
features which are typical of the isotropic model. In the 
ordered phase it is the specific heat behaviour that testify 
to the existence of long- wavelength low-energy isotropic- 
like excitations. On the other hand, when the criti- 
cal region is abandoned in the paramagnetic phase, the 
anisotropy looses its fundamental role and a crossover 
towards the isotropic behaviour is observed, at least as 
far as the correlation length is concerned; such crossover, 
however, has a weaker meaning for larger anisotropy, be- 
ing confined, already for ^ = 0.7, to the high temperature 
region where ^ is the order of the lattice spacing and dif- 
ferences between different models become irrelevant. 

We have compared our theoretical results with the 
neutron scattering experimental data for the staggered 
magnetization, staggered susceptibility and correlation 
length of the real compound Rb2MnF4 and found an ex- 
cellent agreement both for the overall temperature be- 
haviour and for the value of the critical temperature. 

We can hence conclude that the experimentally ob- 
served finite temperature transition in Rb2MnF4 is due to 
an easy-axis anisotropy in the intra-layer exchange inter- 
action and that, despite the small value of the anisotropy, 
the compound shows an Ising-like critical behaviour. 
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